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We study the multiband superfluid and the superfluid (SF) to band insulator (BI) transition 
of strongly interacting fermionic atoms in an optical lattice at a flUing of two fermions per well. 
We present physical arguments to show that a consistent mean held description of this problem is 
obtained by retaining only intraband pairing between the fermions. Using this approach we obtain 
a reasonable account of the experimentally observed critical lattice depth for the SF-BI transition 
and the modulated components of the condensate density, and make predictions for the lattice 
depth dependence of the quasiparticle gap which can be tested in future experiments. We also 
highlight some interesting features unique to cold atom superfluids within this intraband pairing 
approximation; for instance, the pair field is forced to be uniform in space and the Hartree field 
vanishes identically. These arise as a result of the fact that while the pairing interaction is cut off 
at the scale of the Debye frequency in conventional superconductors, or at the lattice scale in tight 
binding model Hamiltonians, such a cutoff is absent for cold Fermi gases. 



I. INTRODUCTION 

Studies of cold atomic gases have led to significant 
experimental [l| and theoretical progress in our un- 
derstanding of the crossover regime between a weakly 
paired BCS superfluid and a strongly paired BEC super- 
fluid. Recent experiments have also studied a strongly 
interacting fermion superfluid (SF) in an optical lattice. 
Upon imposing the optical lattice potential on a SF of ^Li 
atoms with near-unitary scattering, at a density of two 
atoms per lattice well, it was found [sl that the SF re- 
tained its coherence at weak to moderate lattice depths 
as deduced from the interference peaks in the momen- 
tum distribution of Feshbach molecules. Beyond a cer- 
tain critical lattice depth, however, the interference peaks 
disappeared, signalling a possible SF to band insulator 
(BI) quantum phase transition. Motivated by these ex- 
periments, theoretical studies 0, Q have focussed on the 
BI states of strongly interacting fermions in deep lattices, 
and shown that such Bis can be unstable to Cooper pair 
formation below a critical lattice depth. Interestingly, a 
toy model for such an interaction dependent SF-BI tran- 
sition was proposed and studied some years ago [6] in the 
context of the cuprate superconductors. 

In this paper, we focus on the problem of understand- 
ing the multiband SF, and SF-BI transition, of fermions 
in a three dimensional (3D) optical lattice, at a density 
of two fermions per well, using a mean fleld intraband 
pairing approximation which is equivalent to Anderson's 
idea of pairing exact time-reversed eigenstates . This is 
known to be an excellent approximation in the context of 
disordered solid-state superconductors [§] and in metallic 
nanoparticles P|, where pairing between single-particle 
states which are not related by time-reversal symmetry 
gets suppressed by quenched randomness. In the present 
context, we argue that the intraband pairing approxi- 



mation provides an internally self-consistent mean-field 
approach and that interband pairing is expected to be 
dynamically suppressed. Using this theory we obtain sev- 
eral interesting results, summarized below. 

On the theoretical front: (i) We point out that there 
exists a significant discrepancy between the earlier theo- 
retical estimates [1, Q of the critical lattice depth for the 
SF-BI transition, and we provide a physical resolution of 
this discrepancy, (ii) Our mean field approach provides a 
numerically tractable route to understanding properties 
of the SF phase which was not considered in earlier the- 
oretical work 0,111. (iii) Finally, we show that there are 
a number of ways in which the cold atom multiband SF 
is distinct from multiband solid state superconductors. 
These differences arise from the fact that in solid state 
systems there is a cutoff in the pairing interaction, set by 
the Debye frequency in a conventional superconductor or 
by the lattice cutoff in model tight-binding Hamiltoni- 
ans, which is absent in cold atom SFs. More generally, 
this implies that results obtained on such effective lat- 
tice models in the context of cold atom Fermi superfluids 
must be interpreted with caution. 

On the experimental front: (i) As seen from Fig. 1, 
we find that the SF-BI transition, at unitarity, oc- 
curs at a critical lattice depth V^^"^^ ~ ^En-, where 
Efi = Ti'^h^ / {2mab\) is the atomic recoil energy, with 
6l — A/2 being the lattice spacing (A being the laser 
wavelength) and rua being the atom mass. Our result 
is somewhat larger than the experimental estimate Q of 



yicrtt) ^ ZEr, but is reasonably close, given that our 
approach is a mean-field theory, which neglects fluctu- 
ations (but see below). Our result for V^^"^^ coincides 
with the result of Ref. [ij but follows from a set of in- 
ternally consistent assumptions. We also show why the 
result of Ref. while obtained within a valid formalism. 
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leads to a significant overestimate of tfie critical lattice 
depth for the SF-BI transition at mean field level, (ii) 
Fig. 1 shows the uniform and modulated components of 
the condensate density. While the uniform component 
decreases monotonically with the lattice potential, the 
different modulated components are nonmonotonic since 
they are induced by the lattice potential. The maximum 
value of the first order peak occurs at about half the crit- 
ical lattice depth and the peak value is about a factor of 
five smaller than the uniform component, which is in a 
good agreement with the experimental data (iii) As 
seen from Fig. 2, the excitation gap for fcrmionic quasi- 
particles is always nonzero at low temperature, both in 
the SF and in the BI with the minimum excitation gap 
occurring for lattice depths slightly smaller than the crit- 
ical lattice depth. 

The paper is organized as follows. Section [Til contains 
a theoretical discussion where we develop our mean-field 
theory and show that the intraband pairing ansatz gives 
an internally consistent description of the physics at den- 
sity of two atoms per well. Section lllll explores the exper- 
imental consequences of our mean-field theory. Specifi- 
cally, we present results for the critical lattice depth for 
the SF-BI transition, the dependence of both uniform and 
modulated components of the condensate density in the 
SF phase on the lattice depth and compare these results 
with the experiment. We also compute the quasiparticle 
excitation gap as a function of the lattice potential depth 
and find a nonmonotonic dependence which can be tested 
in future experiments. We conclude in Section ITVl with a 
brief summary of the results. 



II. MODEL AND MEAN FIELD THEORY 

Let us consider fermions moving in a 3D lattice poten- 
tial 



Vr 



,2ttX, ,2ttY , ,2ttZ, 
cos(— — ) + cos(— — ) + cos(— — ) 
ol Ol Ol 



(1) 



with lattice constant = A/2, formed by three pairs of 
counterpropagating laser beams with wavelength A. Let 
us use units where we measure momentum (distance) 
in units of and energies in units of the sin- 

gle atom recoil energy, Er = ir'^h? /2mab\. Solving the 
(dimensionless) single-particle Schrodinger equation, we 
find band energies e„k and the Bloch functions u„k(r). 
Expanding the fermion field operator in the Bloch basis 
as ip{v) = I]nk V'«kWnk(r)e''' '', where N is the total 
number of unit cells, and including the interactions, leads 
to the Hamiltonian 

if = XI ^»kV'lk>„k<T - j <f'rij\{v)^\{r)'4j^{r)i>^{r). 

(2) 

Here ^„k — £nk — M (/^ is the chemical potential), 
= 9s/{b\l^Ft) is the (dimensionless) contact inter- 
action which can be tuned with a magnetic field, and 175 



determines the s-wave scattering length, ag, via 



4:7rh as 



55(A) 



(27r)3 n^q^ 



(3) 



As indicated, gs must depend on the ultraviolet cutoff A 
to recover the measured (low energy) scattering length 
as- We ensure that A is large enough for physical prop- 
erties to have converged to a cutoff-independent value in 
our numerical calculations. An important point to note 
is that gs (and hence ws) scales as ~ 1/A when A ^ 00. 

To proceed with the many-body problem, we first 
follow the standard Bogoliubov-de Gennes method [loj 
and introduce mean fields p{r) = X]o-(^<T('')^cr('")) ^^'^ 
A(r) = ■u;s(V'|(r)V'|(i")), where the expectation values 
are evaluated in the mean-field ground state which is to 
be determined self-consistently. Previous work [1, Q ig- 
nored the Hartree mean field p{v); we first show that 
this is justified in the cold atom superfluid despite the 
inhomogeneous density induced by the optical lattice. 
In order to see this, we note that the inhomogeneous 
Hartree shift is simply a renormalization of the optical 
lattice potential as Ve//(r) — VL{r) — wsp{r)/2. In the 
limit that the cutoff A 00, ws ~ 1/A, while the den- 
sity p{r) is certainly finite everywhere. This means that 
Ve//(r) = ^^(r) and that the Hartree terms play no 
role in the limit A 00. We will later see that this 
same argument fails for the pairing potential A(r) since 
('0|(r)'!/'^(r)) is in fact ultraviolet divergent unhke p{r); 
this divergence is precisely compensated by ws ~ 1/A. 

Let us proceed by introducing Fourier modes for the SF 
order parameter modulations as A(r) — J2g ^(g)e*^ "^j as 
well as for the Bloch functions, M,ik(r) = J2g "nk(g)e*^ '", 
where g are reciprocal lattice vectors. w„k(g) can be 
chosen to be real due to time reversal symmetry. We can 
then define matrix elements 



A„„'k = X A(g)F„„/k(g), 



(4) 



where F„„/k(g) = «nk(g + g')Mn'k(g')- The pairing 
matrix element A„„'k corresponds to singlet pairing of 
fermions at momenta (k, — k). This leads to intraband 
pairing for n = n' and interhand pairing for n ^ n' . Both 
types of pairing matrix elements are generally nonzero. A 
mean field decoupling of the interaction term in H leads 
to 



nkcr 

- X [A„„/kV'LT^l'-ki+^™'kV'™'-kiV'«kT -(5) 

nri'k 

The self-consistent mean field state can be obtained iter- 
atively, by diagonalizing Hmf in Eq. ([5]) and using those 
eigenstates to evaluate the pairing field matrix elements 
in Eq. ^ which yields the new Hmf- This is an ex- 
tremely computationally demanding task due to a large 
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number of bands 100) one needs to keep in order 
to converge to cutoff-independent results at unitarity. 
However, if one is interested in the critical lattice depth, 
y{cnt)^ at which the SF to BI transition occurs, the com- 
putations can be significantly simplified by studying the 
linearized Bogoliubov-de Gennes equations, 



A(g) = ^/C(g,g')A(g'), 



(6) 



g' 



where the pairing kernel is given by 



-n«'k(g)-F„„'k(g')- (7) 

Here np is the Fermi distribution function. (Ref.[3] in- 
correctly assumed that this pairing kernel was diagonal 
in g.) One can calculate vjf"*^ by diagonalizing the ma- 
trix 5g_g' — A^(g, g'), as was done in Ref. and finding 
the value of Vl at which its lowest eigenvalue becomes 

negative. This was shown [5j to lead to V^^"^^ ~ 45i?fl 
which is more than an order of magnitude larger than 
the experimental estimate y^'^'"'*-' « iEji. 

The reason for this enormous discrepancy can be un- 
derstood if we look at the real-space structure of A(r), 
corresponding to the lowest eigenvalue of the pairing ker- 
nel A^(g, g') near vjf"^\ by Fourier transforming the cor- 
responding eigenvector. We find that A(r) near vj^'^"*^ is 
very strongly modulated; it is large near the center of each 
well of the optical lattice but is nearly zero in the region 
between adjacent wells. For a filling of two particles per 
well (p = 2), this means that A(r) varies rapidly on the 
scale of the interparticle spacing 1/kp which is 

comparable to the coherence length at unitarity. Am- 
plitude fluctuations of the order parameter will therefore 
be very important for such rapidly varying components 
of A(r) and will tend to suppress such fast modulations. 
Since intcrband pairing matrix elements only arise from 
A(g ^ 0), this also means that the amplitude fluctua- 
tions will predominantly suppress intcrband pairing am- 
plitudes. In short, a mean-field treatment of intcrband 
pairing is internally inconsistent at the atom density of 
two atoms per well: it predicts a significant variation of 
the pairing mean field on short length scales, where am- 
plitude fluctuations are expected to be significant. We 
do expect such a mean field theory to provide a better 
starting point in the case where there is a large number 
of atoms in each well so that amplitude fluctuations on 
the scale of 6^ can be expected to be small. In this case, 
the dominant fluctuations will be quantum phase fluc- 
tuations which can drive a SF-BI transition analogous 
to the superconductor-insulator transition in Josephson 
junction arrays ll|. We will confine our attention here 



to a density of two atoms per unit cell. 

In contrast to intcrband pairing, the intraband pairing 
is stronger when the mean-field potentials are more ho- 
mogeneous, as seen from Eq.Q. Thus one can obtain a 



self-consistent mean-field description of the SF-BI tran- 
sition in our system by just retaining intraband matrix 
elements in Eq.(l5]). This corresponds to Anderson's idea 
of pairing time-reversed eigenstates Q • This is known to 
be an excellent approximation in the context of the or- 
dinary solid-state-based inhomogeneous, in particularly 
disordered, superconductors [sl] and in metallic nanopar- 
ticles [9] . In this case pairing between non-time-reversed 
eigenstates is suppressed by the randomness of the cor- 
responding matrix elements. In our case the suppression 
is dynamical leading, however, to the same final result. 

Under this intraband approximation, the mean-field 
Hamiltonian is obtained by setting A„„/k = A„k^n,n' in 
Eq. (5). Defining the Bogoliubov quasiparticle disper- 
sion as £nk — \/C^k + ^nk' ^^'^ *° mean-field 
equations 



^(g) = ^E^t-h('^)F„„k(g), 



nk 
nk 



2T 



t-nk V 2T 



(8) 



where p = 2 is the atom density per well. The pair- 
ing matrix elements are in turn determined via A„k = 
A(g)_F'„„k(g)- In order to solve these equations we 
need to restrict the number of bands. This is done by 
assuming |g| < 27rA, where A is the cutoff in Eq. ([3]). 



III. RESULTS 

We now turn to some of the key results obtained from 
solving the intraband mean-field equations that we moti- 
vated in the previous section. We begin by pointing out 
some general qualitative ways in which cold atom super- 
fluids differ from ordinary solid state systems although 
the mean field gap and number equations appear to be 
similar for these two systems. We then turn to results ob- 
tained from numerically solving the mean field equations 
for some experimentally relevant observables. 



A. What is special about atomic Fermi superfluids? 

We begin by pointing out an interesting and unex- 
pected consequence of the intraband pairing ansatz — 
namely, the pairing field in real space is completely uni- 
form when we take the limit A ^ oo as we should. This 
can be seen by inspection of the expression for the Fourier 
component of the pairing field at wavevector g in Eq. ([8]). 
As seen from Eq. (3), the coupling constant ws vanishes 
as ~ 1 / A as the cutoff is taken to infinity. This vanish- 
ing of the coupling constant is cancelled by the expansion 
of phase space for pairing, i.e. the number of bands at 
high energy which contribute to the sum in Eq. ([8]), pro- 
vided the matrix element Fnnkis) allows it. However, 
high energy (e„k ^ Vl) states are simply free-particle 
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states for which F„„k(g) ~ ^g,o- This means that only 
A(g = 0) component of the pairing field survives in the 
limit A — * oo. Note that this result is independent of the 
form of the lattice potential and can be expected to hold, 
for example, in the case of the random optical lattice po- 
tential, which can be created by laser speckle [13]. This 
has potentially interesting implications for the disorder- 
driven SF to insulator transition in cold atom SFs. Also 
note that these observations apply only to the cold atom 
SFs; in conventional superconductors, the finite cutoff (at 
the Debye frequency in conventional superconductors, or 
at the lattice cutoff in model tight binding Hamiltonians) 
does allow for spatial variations of A(r) within the ap- 
proximation of pairing time-reversed eigenstates In 
terms of a Landau theory, if we have a uniform superfluid 
order parameter and a periodic density modulation, then 
a symmetry allowed term in the Landau functional such 
as / cfirp{r)\A{r)\'^ will certainly lead to an order param- 
eter modulation with the same period as the density mod- 
ulation. Our results indicate that the coefficient of this 
term must vanish, in the intraband pairing approxima- 
tion, as the cutoff A ^ oo. One final peculiarity of cold 
atom SFs is that the pairing field (and the pairing gap) is 
completely rigid^ i.e. it does not depend on the band in- 
dex n or the momentum k, unlike in multiband solid state 
superconductors. To the extent that interband pairing is 
suppressed, as we have argued is the case for two atoms 
per well, our above observations imply that any attempt 
to model this multiband cold atom SF by imposing a 
finite band cutoff, as is customarily done in deriving ef- 
fective tight-binding lattice models, will generally lead to 
incorrect results. We now turn to results obtained from a 
numerical solution of the intraband mean-field equations 
for observables of experimental interest. 




2Vl/Er 



FIG. 1: Lattice depth dependence of the condensate density 
in the SF phase for different scattering lengths. nc(0) denotes 
the uniform component of the condensate density while nc(l) 
and nc(2) correspond to modulated components at wavevec- 
tors g — [2n, 0, 0] and g = [2tt, 2tt, 0] respectively. 



C. Condensate density 

It is important to realize that uniformity of the pairing 
field does not imply uniformity of the condensate density, 
which is what was measured in the experiment of Ref. 
Q. The condensate density is obtained from the two 



particle density matrix [14l |. In a lattice environment. 



this condensate density has Fourier components 



nk 



(g) (V'ikrV'. 



I-ki) 



(9) 



B. Critical lattice depth for the SF-BI transition 



The sum over the band index n in Eq. ^ is well defined in 
the limit A —^ oo and all Fourier components of ric are in 
general nonzero. As shown in Fig. 1, the uniform compo- 
nent nc(g = 0) decreases monotonically with increasing 
lattice depth, while the modulated components have a 
nonmonotonic dependence on Vl as observed in the ex- 
periment. The ratio of the leading modulated component 
to the uniform condensate density is in a good agreement 
with the data. 



We begin by investigating the critical lattice depth for 
the SF-BI transition at atom density p — 2. For fermions 
with unitary scattering, the SF order parameter vanishes 
at V^"*'^ ~ AEr, close to the value given in Ref. [4] ob- 
tained by simply assuming, in effect, that the order pa- 
rameter is uniform. The results and arguments presented 
in this paper justify such an assumption. As we already 
mentioned, experimentally it appears that y^"*''^ "AEr 
(note that Chin et al. 3] quote their results in units of 
the molecular recoil Efi/2). Part of the reason for the 
remaining relatively mild discrepancy with the experi- 
ment might be quantum and thermal fluctuations beyond 
mean-field theory, which are known to be somewhat im- 
portant at unitarity even in the absence of the lattice 
One must also keep in mind, however, that there is 
no good experimental estimate of the temperature after 
the lattice is turned on 



D. Quasiparticle gap 

The SF state as well as the insulating state have a gap 
to quasiparticle excitations. On the SF side, the quasi- 
particles are Bogoliubov quasiparticles while the quasi- 
particles in the insulating phase correspond to the origi- 
nal fermions. As shown in Fig. 2, the minimum excitation 
gap, 2Eg, in the SF phase decreases with Vl (tracking 
2A(0) over a wide range of lattice depths) whereas it in- 
creases with Vl on the insulating side. This leads to a 
minimum in the excitation gap for lattice depths close 
to but below yj^'^"*^ . This nonmonotonic variation of the 
gap, shown in Fig. 2, could be tested in experiment. Note 
that there is a small window of lattice depths, close to 
but below vl^'^"^\ where the quasiparticle gap is slightly 
less than the gap of the underlying BI. This arises be- 
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FIG. 2: Minimum excitation gap, 2Eg, as a function of lattice 
depth in the SF and BI phases at low temperature for scat- 
tering lengths corresponding to unitarity (open circles) and 
for as/bh ~ —1 (open squares). Solid line indicates the in- 
sulating gap of the band problem. This gap corresponds to 
what would be measured in typical experiments probing the 
gap which will create two Bogoliubov quasiparticles in the SF 
or a particle-hole pair in the BI. 

cause of particle-hole symmetry breaking and would not 
occur in particle-hole symmetric toy models 



IV. CONCLUSIONS 

In this paper, we have discussed the multiband SF 
phase of strongly interacting fermions in an optical lat- 
tice. Our results for the critical lattice depth of the SF-BI 
transition and the dependence of the Fourier components 
of the condensate density on the lattice depth are in a 
good agreement with experimental data of Ref . [3] . Our 
prediction for the quasiparticle gap in the SF and insu- 
lating phases could be tested in future experiments. We 
have also obtained new results of general relevance to cold 
atom SFs. In particular, the uniformity of the pairing 
amplitude under time-reversed-eigenstate pairing condi- 
tions is an interesting and surprising result which can 
be expected to have a strong effect, for example, on the 
nature of the disorder-driven SF to insulator transition 
in cold atom SFs. This will be explored in future work. 
Other directions for future research include a study of 
quantum fluctuations in the multiband SF, dynamical 
instabilities of such fermion SFs and the study of 
SF-insulator transitions in other lattice geometries fl^. 
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